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A relativistic description of pion-nucleon scattering based on the four- 
dimensional Bethe-Salpeter equation is presented. The kernel of the equation 
consists of s- and li-channel nucleon and A(1232) pole diagrams, as well as 
p and a exchange in the i-channel. The Bethe-Salpeter equation is solved by 
means of a Wick rotation, and good fits are obtained to the s- and p-wave irN 
phase shifts up to 360 MeV pion laboratory energy. The coupling constants 
determined by the fits are consistent with the commonly accepted values in 
the literature. 



I. INTRODUCTION 

Pion-nucleon (tcN) scattering is an important example of a strong interaction, and as 
such plays a significant role in many nuclear reactions involving pions, the most interesting 
example in recent years being pion photoproduction. It is generally accepted that the fun- 
damental theory of strong interactions is Quantum Chromodynamics (QCD), and therefore 
a theory describing ttN scattering should ideally be derived from QCD. However, due to 
the non-perturbative nature of confinement, QCD has not been amiable to solutions for low 
and intermediate energies, and hence it is necessary to use an effective theory. In principle 
the effective theory should be as close as possible to the fundamental one, and so should 
satisfy the same symmetries, in particular chiral symmetry, which is known to be impor- 
tant for low-energy physics. Therefore, in place of the QCD Lagrangian a chirally-invariant 
hadronic Lagrangian is used, where the degrees of freedom are mesons and baryons rather 
than quarks and gluons. For low energies it is expected that the detailed quark structure 
of hadrons is relatively unimportant, and that it is only at very high energies that explicit 
quark degrees of freedom are essential. The great success of meson-exchange models for 
nucleon-nucleon (NN) scattering in the description of NN phase shifts is an example of 
this. 

A number of dynamical models of nN scattering have been developed over the past few 
years. Most begin with a potential which is iterated in a Lippman-Schwinger-type equation 
to give the scattering amplitude, from which the phase shifts and observables are obtained. 
This method ensures that two-body unitarity is respected, and that multiple scattering 
effects are taken into account. The simplest models use separable potentials fTpH, in which 
the parameters have no physical meaning, and furthermore, different sets of parameters 
are used in each partial wave. While these models provide good descriptions of the ttN 
phase shifts, they provide no information about the interaction process. An alternative is 
to derive a potential from a Lagrangian which describes the couplings between the various 
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mesons and baryons. In tree-level models P-pl, the potential is unitarized using the K- 
matrix approximation, which relies on the assumption that for low energies the .fT-matrix is 
equal to the potential. While this method can provide a good description of the ttN data 
at low energies, in order to cover a larger region of energies and to be able to investigate 
the nature of resonances, the potential must be iterated to all orders. This has been done 
recently JIT -18] in models based on three-dimensional (3-D) approximations to the Bethe- 
Salpeter equation |19||. Phenomenological form factors are required at the vertices in order 



to provide convergence, and so in general more free parameters are needed than in tree- level 
models, which can give good results without any form factors ||. 

There are several other approaches to ttN scattering in the literature which give results of 
a similar quality. The meson-exchange model of Jade [2(J uses the solitary boson exchange 
method to regularize self-energy diagrams instead of using form factors, and gives good 
descriptions of both ttN and NN scattering with the same set of parameters. Sato and Lee 
constructed a meson-exchange model using an effective Hamiltonian PT|, while Ellis and 



Tang used chiral perturbation theory. Quark models have also been used to describe 
ttN scattering f23^p5|. Fuda 0] developed a Poincare invariant front form model of ttN 
scattering, in which the potentials were assumed to be seperable. This model was later 
extended |jj to a pion laboratory kinetic energy of 1 GeV, giving an excellent description of 
the phase shifts and inelasticities in the s-, p- and <i-waves. 

The exact ttN <— tcN amplitude for a given Lagrangian can in principle be obtained 
from the full BS equation, with fully dressed propagators in the nN intermediate state, and 
a potential consisting of all 1- and 2-particle irreducible connected diagrams. By definition 
this exact amplitude would satisfy both crossing and chiral symmetry, and have the correct 
one-body limit (an equation is said to have the one-body limit if it correctly reduces to either 
the Klein-Gordon or Dirac equation if the mass of either of the particles becomes infinitely 
heavy). Unfortunately, at present it is impossible to construct the potential for the full 
BS equation, as it would contain an infinite number of Feynman diagrams. Consequently 
the potential and propagators must be approximated in some way. The simplest and most 
commonly used approximation involves replacing the dressed propagators in the intermediate 
states with bare propagators with poles at the physical masses, and truncating the potential 
so that it contains only the lowest-order diagrams, i.e. the tree-level diagrams. The resulting 
equation is generally referred to as the ladder BS equation. There are some problems with 
this equation, such as it does not have the correct one-body limit [26[]. Also some symmetries 
present in the approximate potential, such as crossing and chiral symmetry, may be violated 
in the solution of the BS equation |27| . 



The BS equation for the nN scattering amplitude is a covariant four- dimensional integral 
equation, in which the integration is over the relative energy and relative momentum of the 
ttN intermediate state. It is the presence of the relative energy as an integration variable 
that is responsible for the complicated singularity structure of the kernel. Frequently the BS 
equation is reduced to a 3-D integral equation in order to avoid the difficulties involved in 
the handling of the singularities of the kernel. This is achieved by approximating the kernel 
in such a way that the integration over the relative energy can be carried out explicitly, 
resulting in a 3-D integral equation. There are an infinite number of 3-D reductions of the 



BS equation [28], all of which satisfy relativistic elastic two-body unitarity. There is no over- 



whelming reason to choose one particular 3-D approximation over any other. Some of the 
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3-D equations are chosen so as to overcome problems caused by approximating the full BS 
equation with the ladder approximation. For example, the "smooth propagator" [29], used 
by Pearce and Jennings JT1J and Jade |[2Q|| , attempts to restore chiral symmetry. However, 
other problems are introduced when the dimensionality is reduced from 4-D to 3-D. Pasca- 
lutsa and Tjon [[H] have shown that when the nucleon self-energy is calculated using most 
of the commonly used 3-D propagators, there are differences in the renormalization between 
the positive and negative energy states, which is an indication that charge conjugation and 
CPT symmetries are violated. Another problem is that there are significant differences in 
the half-off-shell amplitudes as calculated by different 3-D reductions of the BS equation 
even when they all give the same results on-shell. This could have significant implications, 
in that it is the off-shell behaviour that is important when the nN amplitude is used as 
input into the calculation of other nuclear reactions, such as pion-nucleus scattering, pion 
production in nucleon-nucleon collisions, or pion photoproduction. 

In the present investigation we describe a relativistic model of nN scattering in which the 
BS equation is solved directly in four- dimensions. In this way we can avoid the ambiguities 
encountered in reducing the BS equation to three- dimensions. This we hope will give us 
coupling constants that could be compared to those extracted from QCD models, as it has 
been shown that there are considerable differences in coupling constants obtained using 3-D 
equations as compared to the BS equation [30]. More important is the fact that our tcN 



amplitude can be gauged to give a photoproduction amplitude that satisfies unitarity 
and gauge invariance. 

To our knowledge there are currently no models based on the BS equation that give the s- 
and p-wave ttN phase-shifts in good agreement with the empirical data. The only previous 
meson-exchange model for ttN scattering to have used the 4-dimensional BS equation was 
that of Nieland and Tjon Since the potential consisted only of the w-channel nucleon 
pole diagram, this model could only give the P33 phase-shifts in agreement with experiment. 
In order to get a good description of all the s- and p-wave phase shifts, it is necessary to 
include in the potential s- and ^-channel N and A (1232) poles, and in addition t-channel p 
and a exchange diagrams. All of the recent irN meson-exchange models have included these 
diagrams, although Schiitz et al. replaced the p and a exchange diagrams with correlated 
two-pion exchange [ 15 1 . The model of Schiitz was later extended |16] to include the coupling 
to the ir A, rjN and crN channels, and included the iV*(1535) pole diagrams in the potential. 
Pascalutsa and Tjon included the Roper resonance |p7j , and later also included the Sn and 
D ri resonances as elementary particle poles in the potential |T5 |. We do not include the 
Roper or any other higher baryon resonances, or coupling to channels other than ttN, as 
these contributions are expected to be small for elastic ttN scattering and at pion energies 
below the Roper resonance. Here we are mainly interested in ttN scattering below the 
two-pion production threshold (around 360 MeV pion laboratory energy). 

The organization of this paper is as follows. In Sec. |I| we state the Bethe-Salpeter 
equation and give expressions for the dressed vertices and propagators. The Lagrangian and 
form factors used in our model, and the choice of spin-3/2 propagator for the A(1232), are 
discussed in Sec. |TJ. The renormalization procedure is outlined in Sec. [TV|, and we then 
proceed to describe our method of solving the BS equation in Sec. |V|. Here we also discuss 
the analytic structure of the BS equation. In Sec. [\T| we present results of fits to the SM95 
partial wave analysis of Arndt et al. PB[ from threshold up to 360 MeV pion laboratory 
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energy. We compare the coupling constants obtained in the present work to those obtained 
using models of ttN scattering based on three-dimensional reductions of the BS equation. We 
also discuss how the phase shifts are built up from the individual diagrams in the potential, 
and calculate the renormalized pion cutoff mass in order to examine the effect of dressing 
on the nNN form factor. Next we consider a second model which differs from the first in 
the parameterization of the form factors. We conclude this section by presenting results for 
the phase shifts up to 600 MeV. Finally, in Sec. |V11| we present some concluding remarks. 



II. THE BS EQUATION FOR irN SCATTERING 

We consider the scattering process 

*W + Wn) — *M + N{p N ) , (1) 

where p n and p^ represent the incoming 4-momenta of the pion and nucleon, while p' n and 
p' N are the outgoing 4-momenta. The total 4-momentum P is given by 

P = Pn+PN=p' 7T +p' N , (2) 

and the relative 4-momenta in the initial and final states are defined as 

q = \(jpN- Ptt) , q' = Wn - p'n) • (3) 

In the centre-of-mass (CM) frame the total 4-momentum is related to the total energy y/s 
by P = (\/s, 0). In addition, the Mandelstam variables u and t are given in terms of the 
relative momenta as 

u = (q + q')\ t=(q-q') 2 . (4) 



Having defined the kinematics, the Bethe-Salpeter equation [pl| for the irN <— ttN 
amplitude T(q', q; P) can be written as 

T(q', q; P) = V(q', q; P) - J d 4 q" V(q', q"; P)G nN (q"; P)T(q", q; P) , (5) 

where V is the potential. The two-body ttN propagator G^n is the product of the pion and 
nucleon propagators, i.e. 

r (a . P) _ 1 m + 4 + iriN ( ) 

^N[q,f) (p/2 _ q y _ ^ + %e (p/2 + q y _ m 2 N + ie > W 

where m n and are the charge-averaged pion and nucleon masses. In principle both the 
nucleon and pion propagators in the tcN intermediate states should be dressed. However, 
since here we are only requiring that two-body unitarity be maintained, we have replaced 
the dressed nucleon and pion propagators with bare propagators with poles at the physical 
nucleon and pion masses. 

The potential V(q', q; P) is constructed from the sum of s-, t- and li-channel pole diagrams 
(see Sec. |T|). It is well known that when an s-channel pole is included in the potential of the 
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ladder BS equation, or one of its 3-D approximations, the solution of the BS equation also 



contains an s-channel pole diagram, in which the propagator and vertices are dressed p4 
The potential can be divided into the sum of non-pole and pole contributions, 



V{q', q; P) = V NP (g', q; P) + £ T^ B (q>; P) d$\P) T%(q; P) , (7) 

B 



where T^ ^ is the bare irNB vertex, and d B is the bare propagator for baryon B. At 
present the only baryons we include in the potential are the nucleon and A(1232), and so 
we have B = N, A. The pole part of Eq. (|7]) consists of the sum of s-channel baryon 
pole diagrams while Vnp, the non-pole part of the potential, contains the u- and t-channel 
exchange diagrams. With the potential having this form, the solution of the BS equation 
can be written in a similar way as 

TV, q; P) = T NP (g', q; P) + £ Tl NB (q'; P) d B (P) T wNB (q; P) , (8) 

B 

where IV/vb and ds are the dressed nNB vertex and dressed baryon propagator respectively 
The non-pole part of the T-matrix Tnp is the solution of the BS equation with a potential 
consisting of the sum of the u- and t-channel poles, i.e. 

TkV, g; P) = VW(g', ?; P) - j^y 4 J d 4 q" V np (q', 4'\ ^)GW(g"; T)T NP (g", g; P) . (9) 

Equations (^), ([7p, (|8p and ([]) are shown diagrammatically in Fig. |T|. 

The dressed ttNN vertex is given in terms of the bare vertex and non-pole part of the 
T-matrix as (see Fig. ||) 

IW(g; P) = T^ NN (q; P) - * | d A q"Yf NN {q"- P)G 7rN (q"; P)T NP (g", q; P) . (10) 

The bare and dressed nucleon propagators are 

dP(P) = [f-m£?+ze] -1 , (11) 
d N (P) = [f- m W-S Ar (P)+^]" 1 , (12) 

where m$ is the bare nucleon mass. Also, the nucleon self-energy Sat(P) is given by 

- iE N (P) = -J-^Jd A q rgUg; P)G*n(<1\ PH NN (q; P) , (13) 

and is illustrated in Fig. |3|. Making use of Eq. fllPf) , the nucleon self-energy can be written 



as 



-%z N (P) = --^ / ArS w (g;P)G^(g;P)rS t 7V (g;P) 



I 
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(2tt)« 



J d i q'd i qvf NN {q'-p)GMq l \P)TMq\q\P)GMq\P)^ ) NN{q-,P) ■ (14) 
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The first term in the right-hand-side of Eq. ( |i4|) corresponds to the one-pion loop dressing 
of the nucleon, while the second term is the contribution arising from the iteration of the 
non-pole part of the potential in the BS equation (we refer to this term as the non-pole 
contribution to the self-energy). Since in the ttN intermediate states we have approximated 
the dressed nucleon propagator by a bare nucleon propagator with a pole at the physical nu- 
cleon mass, we avoid having to solve a non-linear Schwinger- Dyson equation for the nucleon 
self-energy. 

There are similar expressions to Eqs. (|1(J) and ([13]) for the dressed 7riVA vertex and the 
A self-energy, however we do not calculate them explicitly, for reasons given in Sec. [TV]. 



III. THE POTENTIAL 

The potential, or driving term, which is iterated to all orders in the BS equation to 
obtain an amplitude satisfying two-body unitarity, consists of all the tree-level Feynman 
diagrams contributing to the process nN <— ttN, derived from the interaction Lagrangian 
under consideration. In this section we discuss the Lagrangian used, the coupling constants 
in this Lagrangian, and the choice of form factors which are necessary to obtain convergence. 
We also look at the possible different forms of spin-3/2 propagators that can be used for the 
A. 



A. The Lagrangian 

The tree-level diagrams shown in Fig. [| are obtained from the following interaction 
Lagrangian: 

Ant = 75 7^ T ■ 9^ IpN + 4>A (^V + x Alplv) TlpN " 0" 77 + h.C. 

ATfl jy YYl^ 

+g<rNN^NlpNO- + O O^TT ■ 8^77 , (15) 

where ipN, V'a,; 7r ; Pm> an d 0" are the fields for the nucleon, delta, pion, rho and sigma, 
respectively. The derivative couplings of the pion field to the other mesons and baryons 
ensures that chiral symmetry is satisfied at tree-level. In the 7riVA Lagrangian, T is the 
transition operator between isospin-3/2 and 1/2 states, and xa is a parameter that can be 
adjusted, and its value will be considered when we discuss our choice for the A propagator. 

Ideally, all of the coupling constants should be fixed using information from other sources, 
rather than leaving them as free parameters. Out of all of the coupling constants appearing 



in the Lagrangian given in Eq. (|T5|), the nNN coupling is the best known (although even 
the Ti NN coupling constant is not without controversy |3li)| , |3"E|] ). We use the value advocated 
by the Nijmegen group p7j , i.e. g^ NN /Ait = 13.5. The remaining coupling constants can be 



determined in a variety of different ways, but there are some discrepancies. 

A value for the ttNA coupling constant can be obtained by calculating the width for 
the decay A — > nN. The coupling constant is chosen such that the width is equal to its 
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experimental value. Assuming that the A self-energy is dominated by the one-pion loop 
diagram, it is found that /^jva/^tt = 0.36 gives the correct width. This result makes use of 
the width calculated on the real s-axis, however it has been shown that there are differences 



between widths calculated on the real axis and in the complex plane f38 |. A somewhat 
smaller value for f n NA is found using the quark-model relation p9f 



LnA = 7TF TT— 9-kNN , ( 16 ) 



25 V2m 



N 



which gives /^tva/^ 71 " = 0.21 when the Nijmegen value is used for g n NN- This smaller value 
has also been shown to be consistent with the width of the A J|0| , provided that higher-order 
mesonic corrections to the A self-energy are included along with the one-pion loop term. 



We now consider the p exchange diagram, which depends only on the product g 



9fmTr9pNN, as well as K p , the ratio of the tensor to vector pNN coupling constants. There 
are a number of different ways of determining a value for g p , as originally discussed by 
Sakurai MM. If it is assumed that the isospin-odd ttN scattering length is dominated by 



the tree-level p exchange diagram, then g 2 p /^ = 3.1 is obtained. Alternatively, it can 
be assumed that the p meson couples to both pions and nucleons with the same strength 
(universality), which means that g p = g p7Tn = g P NN- An estimation of g pnn can be obtained 
from the decay p — > 2tt, which gives g 1 jkn = 2.84. A very similar value of g 1 jAw = 2.7 
is given by the Kawarabayashi-Suzuki-Riazuddin-Fayyazuddin (KSRF) relation |42|], which 



is obtained from current algebra and PCAC, and states that m 2 p = 2g 2 7T7T f 2 , where f n = 93 
MeV is the pion decay constant. There is some controversy about the value of k p . The 
vector meson dominance model (VMD) gives k p = 3.7 [ 43| , while Holder and Pietarinen 



found k p = 6.6 from a tttt — NN partial wave analysis p4 

Although the scalar a meson does not seem to exist in nature, t-channel a exchange is 
included in irN models as an effective interaction, representing higher-order processes not 
explicitly included in the potential, such as correlated 27r-exchange in the scalar-isoscalar 
channel. There are therefore no reliable determinations of the mass of the a meson or the 
magnitudes (or signs) of the aNN and cnrn coupling constants. The mass is usually taken 
to be around 4/77,^ to Qm^. 



B. The A propagator 



There is an ambiguity as to the choice of propagator for a particle with spin-3/2, and 
a number of different propagators have been introduced EBTESl. One way of deriving a 



spin-3/2 propagator is to begin with the free Lagrangian for a massive spin-3/2 field [49 



The propagator obtained from this Lagrangian has, in its most general form, a pole part 
and a non-pole part. The pole part is unique while the non-pole part depends on a complex 
parameter A. When the form of the nNA vertex is chosen correctly, the S'-matrix and 
physical quantities are independent of A p(|. Taking A = — 1 gives the simplest form for 



the propagator, and is commonly called the Rarita-Schwinger (RS) propagator: 



$ + m A 
p 2 — m\ + ie 



9pv 



3m z 



(17) 
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It is known that the RS propagator contains background or off-mass-shell spin-1/2 compo- 
nents, along with the spin-3/2 component (see, e.g. Ref. [pl|). This becomes evident when 
the RS propagator is written in terms of spin projection operators, which are denoted by 
Vfj, and are given by |52| , |5T 



(V 3/2 ) aP = g aP - -~/ a lp ~ 3^ (JlaPp + Valpi) , (18) 

(Vll 2 ) a(3 = ^7a7/3 - + ^~2~ (^7aP/3 + Pal&f) , (19) 

(PU\, = P -f, " " (20) 

{Pil\p = (PaP/3 ~ ilccPp) , (21) 

(^12 2 )a/3 = ^/=^ (f'PalfS ~ PaPfs) ■ (22) 

The RS propagator can then be written in terms of the spin projection operators defined 
above, as 

PM = / + Ti- V 3/2 + -a) V% + -L_ (Vlf + V^) . (23) 

The spin-1/2 background could be considered as being unphysical, since the A(1232) is 
known experimentally to be a particle with spin-3/2. 



Williams introduced a propagator proportional to the spin-3/2 projection operator [46 



M = / I ■ V 3 ' 2 ■ (24) 

The Williams and RS propagators are identical when the A is on-mass-shell, since the spin- 
1/2 components in the RS propagator are only present when the A is off-mass-shell. Also, 
when the Williams propagator is used there are no contributions to the irN <— ttN amplitude 
arising from the parts of the 7riVA vertices proportional to xa- 

There have been attempts to fix ia on theoretical grounds. Peccei |43| suggested that the 
choice xa = —1/4 ensures that there is no direct coupling to the spin-1/2 components of the 
RS propagator, but it was later shown |51| that the spin-1/2 components are always present, 
and cannot be removed by choosing a particular value of x&. Nath et al. suggested that 
xa = — 1 should be used if the ttNA vertex is to be consistent with the principles of second 
quantization. 

The 1/p 2 factor in the Williams propagator can cause numerical difficulties |7,13j. Pas- 



calutsa J47[ used the Hamiltonian path-integral formulation to investigate the interacting 
spin-3/2 field and constructed a theory in which there is no coupling to the spin-1/2 com- 
ponents. For the case of the s- and w-channel tree-level amplitudes for nN scattering, the 
nNA vertex corresponds to the usual ttNA vertex with xa = 0, and the A propagator is 
the same as the Williams propagator but multiplied by p 2 /m\. This extra factor of p 2 in the 
numerator fixes the problems caused by the 1/p 2 term in the spin-3/2 projection operator. 

We consider two possibilities for the A propagator: (i) the Rarita-Schwinger propagator 
and general ttNA vertex, with X A clS db free parameter, and (ii) the Pascalutsa propagator 



S 



and vertex. In both cases the A is treated as a stable particle in the potential, and the 
width is generated dynamically, since the bare s-channel A pole diagram is dressed when 
the BS equation is solved. 



C. Form factors 



The integrals in Eqs. (|5|), (g), QTOD, and ( [U] ) are divergent, and so a regularization 
scheme must be implemented in order to obtain finite results. As is commonly done in 
meson-exchange models, we introduce form factors at each interaction vertex. These form 
factors represent the extended structure of the particles involved, and will ensure that all 
integrals are convergent by suppressing contributions at high momenta. Since at present it is 
not possible to calculate the appropriate form factors directly from QCD, phenomenological 
functions are usually chosen as the form factors, which have no connection to the underlying 
quark dynamics. A consequence of this is that free parameters, the cutoff masses that govern 
the range of suppression, are introduced into the model. 

In meson-exchange models of the NN interaction, the form factors depend only on the 
4-momentum squared of the exchange particles. This cannot be done in ttN models because, 
for example, the form factors would provide no convergence for the s-channel pole diagrams. 
We therefore follow Pearce and Jennings and make the assumption that the cutoff 
function associated with each vertex is a product of form factors that depend on the 4- 
momentum squared of each particle present at the vertex. The abc vertex is therefore given 
by 

Tabc = fabc(ql, ql, <fc) Vabc , (25) 

where V a bc is the coupling operator obtained from the interaction Lagrangian, Eq. (|T5|), and 
the associated form factor is of the separable form: 

fabMl ql ql) = fa(ql) Wb) fMl) • (26) 

The 4- momenta squared of the legs of the vertices are denoted by q 2 , q 2 and q 2 . It is 
conventional to choose the normalization such that f(m 2 ) = 1, where m is the mass of the 
corresponding particle. Therefore, at the unphysical point when all three legs of a vertex 
are on-mass- shell, the corresponding product of form factors is equal to one. 

The scalar functions f(q 2 ) can essentially be chosen in an ad hoc manner, since very little 
is known about the off-mass-shell behaviour of the form factors. One possible choice is the 
multipole form factor: 

where A is the cutoff mass, and n is an integer. Different forms of f(q 2 ) have been used in 
the various models of nN scattering. An example is |i~3| , |T4] 



f(j2) ( (A 2 - m 2 ) 2 r 

MO-I (A2 _ m2)2 + (m2 _ ?2)2 J • (28) 
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The main difference between the two form factors fi and fn lies in the analytic structure in 
the complex q plane (note q 2 = q 2 , — q 2 ). The function fi has two poles along the real q 
axis at 



q = ± yy A. 2 + - iej , (29) 
while fn has four poles in the complex q Q plane, located at 



q = ±jA 2 ±i(A 2 -m 2 )+q 2 . (30) 



The form factors used by Pearce and Jennings JTT| have a very similar structure to fn in 



the complex q plane. We note that fn(q 2 ) always has poles in each of the four quadrants of 
the go complex plane, irrespective of the value of A (except when A 2 = m 2 ). We should also 
mention that fi(q 2 ) with n — 1 could be considered as the propagator for a scalar particle 
with mass A, while fn(q 2 ) with n = 1 corresponds to the propagator for a resonance with 
mass A, and a width proportional to (A 2 — m 2 ). It has been shown that when a form factor 
has the form of a propagator for a resonance with a constant width, there is a violation of 
unitarity at all energies, which only becomes evident when a four-dimensional formulation 
is used [p3[ . 

As an example of the problems caused by this violation of unitarity, we can consider 
the dressed nucleon propagator. If the one-pion loop self-energy diagram is calculated using 
form factors similar to fn{Q 2 ) at each vertex, it turns out that the dressed nucleon has a 
width. This is of course unphysical, since the nucleon is well known to be stable particle. 
As a result, form factors with poles in the complex q plane cannot be used to regularize 
loop diagrams in 4-D models. Therefore, here we use form factors of the form f\, in which 
case the unitarity violations do not occur. 



IV. RENORMALIZATION 

As shown in Sec. O, the bare vertices and propagators appearing in the s-channel pole 
diagrams in the potential become dressed when the potential is iterated in the BS equation. 
A renormalization procedure therefore must be carried out in order to fix the bare parameters 
such that the renormalized quantities are equal to their physical values. Since we include 
form factors at all vertices, the bare masses and coupling constants are finite. 



A. The dressed ttNN vertex 

The TtNN vertex renormalization constant Zin is defined in the usual way: the bare 
ttNN vertex differs from the dressed vertex by only a constant, which is Zuv, when sand- 
wiched between Dirac spinors and all external legs are placed on-mass-shell. Therefore the 
vertex renormalization constant is defined by the relation 

u(P)IW( 9 ; P)u(q) = Z^u(P)T^ NN (q; P)u(q) , (31) 

with P 2 = m 2 N , and q 2 taken such that all three legs of the vertex are on-mass-shell. 
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B. The dressed nucleon propagator 



We now turn to the renormalization of the dressed nucleon propagator. Firstly, Lorentz 
invariance requires that the nucleon self-energy can be written as the sum of a vector and 
scalar part, i.e. 

Sjv(P) =fA(s) +m N B(s), (32) 

where A(s) and B(s) are functions of s only. In order to fix the bare nucleon mass, we 
require that the dressed propagator has a pole at the physical nucleon mass, i.e. 

lim d~ N \P) = 0. 



(33) 



This requirement gives the following expression for the bare nucleon mass in terms of the 
functions A and B: 



m$ = mjy fl — A(m 



B( m %] 



(34) 



The residue of the dressed nucleon propagator at the nucleon pole is defined as the nucleon 
wave-function renormalization constant Z 2 n, i-e. 



lim (f-m N )d N {P) = Z 



2N 



which gives pH] 



where 



J 2N 



A'irni 



N) 



A(m 



d A, N 

ds 



2m 2 N 



(A' (m%) + B'(m%) 



B'(m 2 N ) = -B 



ds 



(35) 



(36) 



(37) 



The above expressions enable the renormalization constants to be calculated easily for the 
case of one-pion loop dressing, since it is not hard in this case to write Sat(P) in terms of the 
functions A(s) and B(s). However, if the non-pole contributions to the nucleon self-energy 
are taken into account, then this would not be the simplest method of calculation. We 
introduce the scalar quantity E^*(s), which is defined as the nucleon self-energy sandwiched 
between Dirac spinors, i.e. 



u(P)Ejv(P)tt(P) 
y/sA(s) + m N B(s) 



since P = in the CM system. Taking note of Eqs. 



and 



constants m N and Z 2 n can be written in terms of S^(s) as 



''V 



[m 



N) 



J 2N 



1 - 2m N —Z™(s) 



-i 



(38) 

, the renormalization 

(39) 
(40) 



In this way we can include both the one-loop and the non-pole contributions to the mass shift 
and wavefunction renormalization Z 2 n in a way consistent with the scattering formulation 
of the BS equation. 
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C. Renormalization of the BS equation 



In the 7rN amplitude there should be factors of \JZ^n for each of the external nucleon legs, 
which result from the application of LSZ reduction |54]] on the ttN <— ttN Green's function. 



However, since we are not performing any explicit LSZ reduction, and we are assuming 
that the nucleon propagators in the kernel of the BS equation are bare propagators with 
physical masses, the factors of y/Zm are not generated, and the ttNN and 7riVA couplings 
are effective coupling constants that should be set equal to the physical coupling constants. 
The only exception is with the s-channel pole diagrams. In this case, the solution of the 
BS equation generates new s-channel pole amplitudes in which the baryon propagators and 
7iNB vertices are dressed. This implies that bare masses and coupling constants should be 
used in the s-channel pole diagrams in the potential. 

In order that the ttN amplitude has a pole at the physical nucleon mass, and that the 
residue at this pole is equal to the square of the physical nNN coupling constant, the bare 
nucleon mass must be fixed by Eq. and the bare ttNN coupling constant fixed using 



9nNN - Z 1N\J Z 2N9nNN • 

Here g^ N is the bare ttNN coupling constant, and g^ NN is the renormalized coupling 
constant, which is set equal to the "experimental" nNN coupling constant by fixing the 
value of g^N N correctly. 

In principle, a similar renormalization procedure should be carried out for the A. How- 
ever, the pole in the T-matrix corresponding to the dressed A occurs in the complex s-plane, 
since the dressed A has a width. Therefore, in order to fix the bare 7riVA coupling constant, 



it would be necessary to analytically continue the BS equation into the complex s-plane [p8 
Rather than doing this, here both and /^a are treated as free parameters. Since the 
P33 partial wave is dominated by the s-channel A pole diagram, the bare A parameters are 
essentially fixed by the P33 phase-shifts. The position at which the phase shifts go through 
90° is determined by the bare A mass, and the width of the resonance is related to the bare 
7riVA coupling constant. 



V. SOLVING THE BS EQUATION 

To calculate quantities such as phase shifts and scattering lengths, we solve Eq. (||), i.e. 
the BS equation with the potential consisting the s- and w-channel N and A poles as well as 
t-channel p and a exchange. In addition, the dressed nNN vertex and nucleon self-energy 
need to be calculated so that the nucleon renormalization procedure can be performed. 



A. Partial wave expansion 

The nucleon propagator present in the irN intermediate states can be separated into 
positive and negative energy components [55], which allows us to write G^n m terms of 
projection operators: 



GWGZ; P) = G au (q , q- s)A + (q) - G™(q , q; s)A"(-q) 



(42) 
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where G uu and G vv are given by 

G uu (q , ^ s ) = ^ ^ /2 + qo -E q + i e + ' (43) 

G* v (q , q;s) = T ?f 1 L , (44) 

£ 9 Vs/2 + g + £<? - «e (V s / 2 - go) 2 - oj 1 + le 



with S g = y q 2 + wi/v an d = \/q 2 + m2 r The positive and negative energy projection 
operators can be written in terms of Dirac spinors as 

A+ (q) = E M r(q)^(q)' ( 45 ) 

r 

A-(q) = -£ u r(q)«r(q). (46) 

r 

The normalization of these spinors is defined by ttj(q)tt r (q) = %f(q)tv(q) = 1- 

The expansion given above in Eq. ((4^) is substituted into the BS equation. Multiplying 
the BS equation from the left and right by Dirac spinors yields a pair of coupled integral 
equations. If we introduce the notation (suppressing the Dirac indices) 

T°V, q; P) = u(q')T(q', q; P)u(q) , (47) 
T™(q', q- P) = v(-q')T(q', q- P)u(q) , (48) 

and similarly for the potential, the BS equation can be written as two coupled equations for 
T Uu and T™: 

T € > u (qi ) q ] P) = V € > u (q>,q-P) 

-T^y E Jd 4 q"V* w "(q', q ";P)G*" w "( q ";P)T*" u ( q ", q ;P), (49) 

V / w"=u,v 

with w = u, v. There is a similar set of two coupled equations for the amplitudes T uv and 
T vv , which are required along with T uu and T vu in the calculation of and Z 2 n- 

We now reduce the number of dimensions from four to two, by removing the angular 
dependence using a partial wave expansion. Including spinor indices again, we can write 
each amplitude in the form 

Atx(q', q; P) = Xl A™' w (q', q; P) X x , (50) 

where w and w' are either u or v, and xx 1S a Pauli spinor. The amplitude A w ' w can be 
expanded in terms of partial wave amplitudes A^iJ as 

A*' w (q' , q'; go, q; s) = N(q', q) E ^ m (q)Af/,J(g , q'; g , g; s)y} /jm (q') , (51) 

it'jm 

where g = |q| , q' = |q'|, and N(q',q) = — (27r) 4 (g'g) _1 . Also, the generalized Legendre 
polynomials are given by 

^jm(q) = E {^m t \m s \ j m)Y irne (q)Xm 3 , (52) 

m,£m, s 
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which are eigenstates of the magnitude of the total angular-momentum operator, J 2 , its 
^-component, J z , the magnitude of the orbital angular- momentum operator, L 2 , and the 
magnitude of the spin operator, S 2 . The partial wave amplitude can be written in terms of 
the original amplitude as 

Affttfo, g'; g , g; s) = j^—^ J dq'dq4 m (q')i^(g^ q'; go, q; s)y iljm (q) . (53) 

Applying the partial wave decomposition to the BS equation and making use of the 
orthogonality of the generalized Legendre polynomials, we obtain 



T|/(g ,g';go,g;s) = ^(g ,g';g ,g; S ) 

/oo roc 
dq'o / dg" VffjM, 4\ C €\ w ' (g ', g"; s)T^(C g"; g , g; s) , (54) 

TetjiWoi 9o, g; «) = v em{^ ?o, g; s ) 

/OO /'OO 
dg ' / dg" Vg^, g'; gff, q"; s)G*" w "(q>>, q"; s)T^(q^ q"; g , g; S ) . (55) 
-oo JO 

w=u,v 



Note that the amplitude Ajpj is diagonal in £, i.e. £ = £', due to parity conservation, 
however amplitudes such as Affij, which involve transitions between positive and negative 
energy nucleon states, are not diagonal in £. Each partial wave amplitude T uu is labelled by 
the orbital angular momentum £, the total angular momentum j, where j — £ ± |, and the 
total isospin /, with / = I or / = |. 



B. Analytic structure 

At this stage it is necessary to examine the analytic structure of the partial wave BS 
equations. We will carry out a Wick rotation []5f|, and analytically continue the ttN am- 
plitude in the q' Q and q' ' variables from the real axis to the imaginary axis. Before doing 
this, we must examine the singularity structure of the kernels of Eqs. ([54]) and (|55| ) in the 
q' ' plane, to make sure that there are no poles or cuts that could interfere with the Wick 
rotation. The residues of any poles present in the first and third quadrants need to be picked 
up, since we rotate the q' ' integration contour from the real axis to the imaginary axis in 
an anti-clockwise direction. The presence of form factors in the potential ensures that the 
kernel is well behaved asymptotically, and as a result there is no contribution from the con- 
tour at infinity. There are three sources of analytic structure that we need to examine: (i) 
the ttN intermediate state, (ii) the potential, in which there are poles from the both the 
exchange particle propagators and the form factors, and (iii) the ttN T-matrix itself. 

The poles of the ttN two-body propagator G uu (qQ,q";s) in the complex q' ' plane are 
located at 



go = ^4(g") = -Vs/2 + Jq" 2 + m 2 N -ie ) (56) 



ql = u;±(g") = T [ ^ q" 2 + m%- * e j , (57) 

corresponding to the positive energy nucleon pole, and the positive and negative energy pion 
poles respectively. G vv (qQ, q"; s) has poles at q' ' = u^, and at 
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q'o = "n(q") = -V^/2 - V<7" 2 +m% + ie, (58) 

which corresponds to the negative energy nucleon pole. If y/s > 2m n it is possible for uo^ to 
be in the first quadrant for < q" < q'^ ax , where g^ ax = Js/A — m\. The positive energy 
nucleon pole can move into the third quadrant for \/s > 2m ^. Therefore, if we stay 
below CM energies of 2m^, we only need to pick up the residues from the positive energy 
pion propagator pole. The residues of G uu (q'o, q"; s) and G vv (q'o, q"; s) at q$ = uj+(q") are 

G ^ q "> ^ = - ,F N , ( R F 1 ' ^ 

IhjqiiLOqii (y/S — tj q ii — 0J q i> + 16) 

GZW; s) = - ™ N . r .„ r • (60) 

We now consider the singularities of all the diagrams present in the potential. Firstly, 
note that the partial wave potentials have the form 

V§M,q';Qo,q"^) = lf) f/x U(x)fr(s,t,u) + ^-P e±1 ( x )f^( s ,t,u)\ , (61) 



N(q', q" 

V$M' 4\ <A, <f\ *) = N( * Qll) f_ x dx ( ^Pt{x)fF(8, t, u) + ^Pe(z)fl°(8, f , u) ) , (62) 



-q^q" 

"1 / n" n' 



where x — q' ■ q" and e q = E q + rriN. The forms of the partial wave potentials corresponding 
to V vv and V vu are very similar to Eqs. (|61| ) and ( j62|) respectively. For the s-channel baryon 
pole diagrams we have f™ w oc l/(s — m%) with i — 1, 2. Therefore, the only analytic 
structure in the q'o complex plane produced by the s-channel pole diagrams is due to the 
form factors on the external pion and nucleon legs, which we will look at shortly. However, 
for the u- and t-channel diagrams we have ff w oc l/(z — in 2 ), where z = u or t, and so the 
functions ff w in this case depend on q' Q \ After carrying out the x integration in Eqs. fl6~I|) 
and (|6"2"D, the partial wave potentials corresponding to the u- and t-channel pole diagrams 
will involve terms such as 

/ (q' + r]q^ 2 - (q' + r]q") 2 - m 2 + ie \ 
\ (q'o + VQo) 2 ~ (q' ~~ VQ") 2 — rn 2 + ie J ' 

where m is the mass of the particle, and rj — 1 (—1) for the ^-channel (t-channel) poles. 
Terms such as these generate logarithmic branch cuts in the q' Q ' plane. For the w-channel 
pole diagrams, the branch points are at 



q'o' = -q' ± (y (?' ± q") 2 + rn 2 B~ iej , (64) 

where B = N or A. The analytic structure for the t-channel exchange diagrams is very 
similar, namely there are branch points at 



qZ = q' ± (xjiqt ±q»y + m\ - ie) , (65) 

where A = p or a. 
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Notice that in Eqs. (0) and ( |65D the positions of the branch points in the q^ plane depend 
on the external variable q' Q , unlike the poles of the ttN intermediate state. Singularities that 
depend on q' are removed by the Wick rotation, since as well as rotating the q^ integration 
contour from the real to the imaginary axis, we analytically continue the ttN amplitude 
to the imaginary axis in q' Q . The branch cuts move away from the integration contour as 
we rotate from the real to the imaginary axis, and in fact the branch points always stay a 
distance m away from the integration contour, where m is the mass of the exchange particle. 

Next we consider the singularities of the form factors. The form factors corresponding 
to the external pion and nucleon legs have n-th order poles at 

<£ = -vV2 ± (vV 2 + A^ - ie) , (66) 
ql = V~s/2 ± ( vV 2 + A 2 - ie) , (67) 

when the functional form given in Eq. fl2"T|) is used. The order of the poles depends on the 
choice for the form factor powers, i.e. n n and n^. To guarantee that these form factor poles 
do not interfere with the Wick rotation, we require that they do not move into the first or 
third quadrants of the q^ plane. This means that we must have \fs < 2A 7T and \fs < 2A^. 
Note that if form factors of the type given in Eq. (p8|) were being used, there would be poles 
from each form factor in both the first and third quadrants. 

As a next step we look at the singularities in the form factors corresponding to the 
exchange particles. For the u-channel exchange diagrams there are branch points at 

ql = -q' ± (y/(q'±q»)* + A% - ie) , (68) 

where again B = N or A. Finally for the t-channel diagrams there are branch points at 

q'l = q' ± (y/tf ± <f ) 2 + Ai - ie) , (69) 

where A = p or a. All of these branch points are removed by the Wick rotation, due to the 
q' Q dependence. 

The remaining source of singularities which must be considered is the T-matrix, i.e. the 
solution of the BS equation. This can be done by looking at what happens as the potential 
is iterated in the BS equation. Singularities in the q' Q plane of T(q', q; P) are generated by 
pairs of poles pinching the q$ integration contour. This analytic structure therefore also 
occurs in the q$ plane of T(q", q; P) which appears in the integrand of the BS equation. As 
the BS equation is iterated a hierarchy of branch cuts are generated. Higher order branch 
cuts arise from the pinching between lower order branch cuts and the singularities of the ttN 
intermediate state and potential. The positions of most of these cuts depends on q' Q , and so 
do not cross the integration contour when both the q' and q^ axes are rotated. However, 
some of the higher order branch cuts do not depend on q' and so are not removed by a Wick 
rotation. If any pair of these cuts in the q^ plane protrude into both the first and third 
quadrants simultaneously, a simple Wick rotation becomes no longer possible. 

Also, singularities pinching the q$ integration contour can produce cuts in the y/s plane, 
i.e. thresholds. If these thresholds are generated by singularities other than those from the 
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form factors, they correspond to physical processes. The lowest-energy physical thresholds 
are at: 

\fs~ = m N + , 

\fs = m N + 2m v , 
\/s = m A + 2m v , 
\/s = m N + m 7T + m a . 

However, if a cut in the i/i-plane is generated by the pinching of a form factor pole and 
another singularity, the threshold is unphysical, since the cutoff mass does not correspond 
to the mass of a physical particle. Some examples of the unphysical thresholds generated 
by the BS equation include: 

yfs = A N + 77V , 

\/s = m N + A T , 
y/s = A N + A n , 
yfs = m N + m w + A a . 

We therefore have to make sure that the cutoff masses are chosen to be large enough so that 
the unphysical thresholds occur above the highest CM energy for which we will solve the 
BS equation. 

In summary, the BS equation can be solved for ttN scattering with our choice of form 
factors using a Wick rotation, provided that the cutoff masses are not too small. There are 
three conditions on the minimum values of the cutoff masses: (i) there are no form factor 
poles in the first or third quadrants of the q$ plane, (ii) any cuts produced by the pinching 
between form factor poles and other singularities are also not in the first or third quadrants, 
and finally (iii) all unphysical thresholds, which are generated by the form factors, are far 
away from the energy region in which we are interested. The minimum values for the cutoff 
masses that can be used in the BS equation are given in Table |. Here we have assumed 
that the two-pion production threshold is the maximum CM energy for which we will solve 
the BS equation. 

If form factors of the type given in Eq. (|28| ) are used, it is not possible to prevent the 
form factor singularities from interfering with the Wick rotation by making any appropriate 
choices for the cutoff masses. There will always be poles from the form factors in each 
quadrant of the plane. 

C. Wick rotation 

Having looked at the analytic structure of the kernel of the partial wave BS equation, 
we are now in a position to perform a Wick rotation on Eqs. fl54|) and (0), by making the 
substitutions 

q'o - »V , <Zo -> < , (70) 

and picking up the residues from any poles in the first or third quadrants of the complex 
plane. After Wick rotation, the partial wave BS equation becomes a system of four coupled 
integral equations. The first two equations for the half-off-shell T-matrix are 
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Tffiiiq'o, q'; go, q\ s) = V$(iq r , q'; g , q; s) 

/oo roc _ _ 

dq'o / dq" 4\ <, q"\ s)G iD w ' «, g"; s)Tf'^{iql q"; g , q; s) 

-oo «/0 

W"=U,V 

- E / <8X dq"V^(iq' , g'; u+(q"), q"; s)G?^ w "(q"; s)Tf^(uj+(q"), g"; g , g; s) , (71) 

Te^ji^q'oi q'i 9o, = Vefjiik'o, q; qo, q; s) 

/OO /" OO _ _ 

/ d?" V^(*g , <?'; 9"; w ' («& g"; s)Tf${iql q"- g , g; a) 
-oo </0 

W"=U,V 

- E (72) 

w"=u,v 

In the above we have put the pion and nucleon in the initial state on-mass-shell. The on-shell 
relative momenta are denoted by go and g, and are given by 



io = 2 ( v q2 + Wj2 n ~ v<? 2 + m l) > ( 73 ) 



and 



/[s- (mjy -hm^) 2 ]^ - (m N - m n ) 2 ] 
9"V S • (74) 

There are two additional equations (usually referred to as the "auxiliary equations" ) which 
are necessary in order to have a closed system of equations to solve: they are Eqs. (|7l| ) and 
(fT2|) with ig replaced with u;+(g'). 

Finally, it is necessary to look at each term present in the potentials to check whether 
there are any remaining singularities after Wick rotation. For energies above the pion pro- 
duction threshold, the w-channel nucleon pole present in the potential in the one-dimensional 
parts of the auxiliary equations develops an imaginary part. This is due to a logarithmic 
singularity moving into the integration region, and must be handled carefully to ensure 
numerically stable results. We do this by carrying out a subtraction similar to Ref. |57| . 



There are no additional singularities caused by the form factors below the two-pion produc- 
tion threshold, provided the cutoff masses are chosen to be larger than the values given in 
Table | 

Above y/s = 2(m -+m 7r ) a cut in the q$ plane, generated by the pinching of the integration 
contour between the positive-energy pion pole and the a meson propagator, moves into the 
1st quadrant. Therefore, above this value of the CM energy, it would become necessary to 
take this additional singularity into account when carrying out the Wick rotation. Here we 
consider CM energies below yfs = 2(m CT + m n ). 



D. Calculation of the phase shifts 

In order to determine the nN phase shifts, the on-shell amplitude T^j(q , q\ q Q , q; s) needs 
to be calculated. This is done by analytic continuation of the half-off- shell amplitude to the 
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on-shell point. In practice, the on-shell T-matrix is obtained using Eq. (ffTD, with both the 
incoming and outgoing particles on-mass-shell. 

The ttN two-body propagator G" e "(g"; s) has a pole when the pion and nucleon are 
propagating on-shell, which occurs for q" = q. This pole is related to the two-body unitarity 
cut. For energies above the 7rN threshold we need to take care of this pole so as to obtain 
equations that can be solved numerically. We achieve this by writing G^ s (q"; s) in terms of 
a principal-value part and an imaginary on-shell contribution, i.e. 

GSfa I ») - - 2Ea „^„ yi _ Eq „ _ —) + ^ Hi - 1) , (75) 

where V denotes that the principal-value prescription should be used when the q" integration 
is performed. 

The ttN phase shifts and the inelasticities r\iji are obtained from the on-shell partial 
wave T-matrix using 



T ejI {q , q; go, q; s) = j — . 76 

m^ir \ 2iq 



The behaviour at threshold is more conveniently described in terms of scattering lengths 
and volumes, which are defined by the effective range expansion: 

q 2e+1 cot S ejI = — + I r m f + . . . . (77) 

Ct£jl I 

Here a^i is the scattering length and r^j is the effective range. 



VI. NUMERICAL RESULTS 

A. Fits to the empirical ttN data 

We begin this section by listing the free parameters in our model. The five cutoff masses 
(Ajv, Aa, A,,., A p and A CT ) are free, but are constrained to be larger than the minimum values 
given in Table |. Due to the uncertainty in the values of the coupling constants other than 
9wNN, we permit f n NA, xa, 9 p , k p and ga-mrQuNN to vary freely. The mass of the a meson, the 
bare A mass and the bare tcNA coupling constant /^a are a ^ so ahowed to vary freely, 
although and are essentially fixed by the P 33 phase shifts. The bare nucleon mass 
rripf and the bare nNN coupling constant g^N are n °t f ree parameters, but are determined 



by the renormalization procedure outlined in Section [TV|. We use the value n a — 1 for all 
the form factor powers, however it turns out that this choice is not crucial to the quality of 
the fit. 

The free parameters are determined in x 2 fits to the s- and p-wave single-energy phase 
shifts up to 360 MeV pion laboratory energy, as well as the scattering lengths and volumes, 
from the VPI SM95 partial wave analysis |33]]. We carry out one fit using the Rarita- 



Schwinger A propagator, and another using the Pascalutsa propagator. The coupling con- 
stants and particle masses for both fits are listed in Table ||. The resulting scattering lengths 
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and volumes are given in Table |T|, and the phase shifts are shown in Fig. |5|. We see that 
the BS equation gives a good description of the irN phase shifts. Notice that the results for 
the phase shifts in the Pn partial wave are better when the Rarita-Schwinger A propagator 
is used. 

As can be seen in Table [II], the values of g a -wK9aNN determined from both fits have negative 
signs, which means that the a contribution is repulsive in the s-waves and attractive in the 
p- waves, as was also found in Refs. [0JT7]] • Note that we have used a low value for the nNN 



coupling constant, i.e. g 2 NN /4ir = 13.5. We have repeated the fits using g\ NN j^n = 14.3 and 
found that the results are very similar to those shown in Fig. |5|. Also, there are no significant 
differences between the coupling constants obtained from the fits using g 2 NN /kK = 13.5 and 
//lv.v/1" 11.3. 

In Table [TV] we compare our coupling constants to those extracted from nN models 
based on the i^-matrix approximation and 3-D reductions of the BS equation. Our coupling 
constants are in general consistent with those obtained from other models of ttN scattering. 
We note that fl N ^/^ is in the range 0.35 to 0.43 for all equations, except for the BS equation 
when the Pascalutsa A propagator is used, in which case the 7riVA coupling constant is about 
twice as large as the commonly accepted value. Our value of xa is similar to those found 
in other models, and they are all in the range —0.41 < ia < —0.11. All the models listed 
give similar values for g 2 , and the values of g 2 all lie in the range 2.5 < gtj^K < 3.36. Our 
values are consistent with the KSRF relation [42] and the value found from the width of the 
pionic decay of the p meson. There is a large range of values for k p , which vary between 
1.44 and 6.6. Our value of k p = 2.66 using the RS propagator is smaller than the VMD 
result of k p = 3.7, and we get a value slightly larger than the VMD result (k p = 4.11) when 
we use the Pascalutsa A propagator. This suggests that within the uncertainty from the A 
propagator, we are consistent with vector meson dominance. 

While we are to get a better fit to the irN phase shifts using the Rarita-Schwinger A 
propagator than when using the Pascalutsa propagator, this does not necessarily suggest 
that the Rarita-Schwinger propagator is the correct spin-3/2 propagator. Other processes 
not included in the present model can give attractive contributions to the Pn partial wave, 
such as the coupling to the 7rA channel or the inclusion of the iV*(1440) resonance into the 
potential. More important is the observation that the different choices for the A propagator 
give rise to differences in the coupling constants. This highlights the importance of having 
a better understanding of how to construct propagators for higher-spin particles. 

Having compared results using two different choices for the A propagator, hereafter we 
restrict ourselves to the Rarita-Schwinger propagator. 



B. Contributions to the phase shifts 

In Fig. ^] we show how the total phase shifts are built up from the contributions of the 
individual Feynman diagrams in the potential. The w-channel nucleon pole is strongest in 
the and P33 partial waves, but also gives important contributions to Pi 3 and P31. The s- 
channel nucleon pole generates the repulsion in the Pn phase shifts, and gives a very small 
contribution to Sn- The w-channel A diagram plays a very important role in all partial 
waves except S31 and P 33 , as it gives a large repulsive contribution to Sn, and gives strong 
attractive contributions to Pn, P13 and P 31 . The s-channel A pole diagram dominates the 
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P33 phase shifts, but also gives an important contribution to the S31 partial wave, and a tiny 
contribution to P31. The contributions from the s-channel A pole to S31 and P31 result from 
the spin- 1/2 components of the Rarita-Schwinger propagator. The p exchange contributions 
are largest in the s- waves, but are also significant in the p- waves except for P 33 . Likewise, a 
exchange is stronger in the s-waves than the p-waves, although the a exchange contributions 
are quite small in all partial waves. 

The Six phase shifts are dominated by p exchange and the ^-channel A pole. By itself, 
the attraction generated by p exchange is far too strong, but is partially cancelled by the 
repulsive w-channel A pole. Except for the w-channel nucleon pole, all of the diagrams 
contributing to the S31 phase shifts are repulsive. The largest contributions to this repulsion 
come from p exchange and the spin-1/2 components of the s-channel A pole. 

The s-channel nucleon pole causes the Pn phase shifts to be negative at low energies. 
The attraction that causes the phase shifts to change sign is dominated by the ^-channel A 
pole and the p exchange diagram. 

The w-channel nucleon pole is strong and repulsive in the P13 and P31 partial waves. In 
fact, by itself, the w-channel nucleon pole almost gives the correct P31 phase shifts. However, 
the p exchange diagram gives a large repulsive contribution to P31, causing the phase shifts 
to deviate strongly from the phase shift analysis. The attraction provided by the w-channel 
A pole to P 31 almost cancels the p contribution completely. Similarly, the P13 phase shifts 
are far too repulsive without the strong attraction produced by the w-channel A pole. 

The P33 phase shifts are of course dominated by the s-channel A pole, with the back- 
ground contribution primarily coming from the w-channel nucleon pole. 



C. Dressing of the ttNN vertex 

As can be seen in Table ||, the cutoff masses we obtain turn out to be quite large. This 
results in the dressing being very significant, as is evident from the large size of the bare N 
and A masses. In view of the significance of the dressing, it is interesting to examine the 
effect of dressing on the irNN form factor. When both nucleons in the bare nNN vertex 
are placed on-mass-shell, the bare ttNN vertex only involves the pion form factor, and is 
given by 

r£(?o, q; s) = fM)K U N N (Qo, q; *) , (78) 

where s = m%, and q\ is the 4-momentum squared of the pion. The bare pion form factor 
is 

Mi) - (t^f )" ■ m 

and we take n n — 1. The renormalized pion form factor can then be introduced as 

(80) 

Pn 

which has the property that f^(m%) = 1, and where q and q are related to the pion 
4-momentum squared by 



TT NN (q ,q-s) 
Vr NN (q ,q;s) 
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2m N 




1 






n 




[™ 2 N 



(81) 
(82) 



In Eq. (|80D we have taken the Pn partial wave of the quantities T^ N (q , q; s) and 
V"]^iv(90)q; s). By comparing the slope of the function f^(ql) at ql = with a monopole 
form factor with cutoff mass A^, we can obtain a value for the renormalized pion cutoff 
mass. We find that = 1.22 GeV, which is softer than the bare pion form factor (recall 
that A,r = 1.77 GeV). This is consistent with previous calculations p^-pT|, which have found 
that dressed form factors are softer than the corresponding bare form factors. 

We can introduce the quantity A,r as a measure of the variation between the renormalized 



pion form factor at q^ 



m; 



and q^ = 0, i.e. 



A„ 



i - /*(o) . 



(83) 



We find that A^ = 1.3%, which indicates that our dressed ttNN form factor is a very slowly- 
varying function of the pion mass. Our value of A^ is somewhat smaller than the value of 
3% obtained using other methods p2 . 



D. Different choice of form factors 

With the choice of form factors we have used so far (hereafter referred to as model I), 
the effect of dressing is significant. We now consider the case where there is a form factor 
only on the pion (referred to as model II), which is arrived at from the parameterization of 
form factors used in model I by taking the limit —>■ oo for h = N, A, p, and o. With 
this choice of form factor there is only one cutoff mass, rather than five, and so the number 
of free parameters is reduced by four. All intermediate states contain the pion propagator, 
and therefore a cutoff function still appears in all loop diagrams to provide convergence. In 
model II the pion form factor is used to vary the off-mass- shell behaviour of the pion. This 
in principle could be constrained by the soft-pion theorems. 

In Table |V| we show three sets of parameters obtained from fits to the on-shell ttN data, 
corresponding to the choices = 2, n n = 4, and n n = 10. The resulting phase shifts are all 
of the same quality as the results using the Rarita-Schwinger A propagator shown in Fig. |5|. 
The coupling constants resulting from the three model II fits are similar to the coupling 
constants obtained using model I, although k p is smaller in model II than in model I, and 
QtrmrQa n n has a positive sign, whereas in model I it is negative. As with model I, there are 
no significant changes to the quality of the fits or values of the coupling constants obtained 
when the fits are repeated using g1 NN /^ = 14.3 as the physical ttNN coupling constant. 

The main difference between the two models is that in model II the effect of dressing 
is not as significant as in model I. The bare baryon masses are much closer to the physical 
masses. The renormalized pion cutoff masses and values of A^ for model II are given in 
Table |V]]. We see that the values of A^ are smaller than in model I, and are close to the 
value of A^ w 0.8 GeV advocated by some authors |p3]] . The values of A n are consistent 
with previous calculations |62j of the difference between f^(m%) and f^(0). 
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E. Above the 2ir production threshold 

In order to see what happens at energies above the 27r production threshold, we show 
the phase shifts up to 600 MeV pion laboratory energy in Fig. |7|. The results from model 
II were obtained using n n = 4. In the Pn and P33 partial waves both models give almost 
identical results, but there are some differences between models I and II in the higher energy 
region in the other partial waves. 

Both the P13 and P31 phase shifts are quite good over the full range of energies, which is 
a reflection of the fact that below 600 MeV there are no resonances in these partial waves 
and the inelasticity is negligible. We require a little more attraction at the higher energies in 
the 531 and P33 partial waves, while a significant amount of additional attraction is required 
for 5*11 above 300 MeV, and for Pq above around 450 MeV. 

It is not unexpected that there are some discrepancies in the higher energy region. Here 
the Sn and Pq partial waves exhibit resonance behaviour not included in the present model. 
There are a number of modifications that could be made to our model in order to improve 
the agreement with experiment for this larger energy range. Firstly, it may be necessary for 
three-body unitarity to be satisfied. Amongst other things, this will involve replacing the 
nucleon propagator in the ttN intermediate states with a dressed propagator. Extending 
the model to include the coupling to inelastic channels and the possible addition of explicit 
nucleon resonances into the potential will be essential at energies above 360 MeV. To fully 
understand the 7rA amplitude at these energies, the coupling to inelastic channels must first 
be included, and if the fit to the phase shifts is still unsatisfactory, explicit bare A* baryon 
poles may need to be included in the potential. 

The 5*ii partial wave would be improved by the inclusion of the coupling to the rjN 
channel, and also Su resonances such as the A* (1535) and A* (1650) may need to be included 
in the potential. The coupling to the ttA and aN channels and possibly the inclusion of the 
A* (1440) resonance into the potential would be necessary to improve the Pn phase shifts 
above 450 MeV. 

VII. CONCLUSION 

In this work we have presented a description of pion-nucleon scattering based on the 
four-dimensional Bethe-Salpeter equation. The kernel of the equation is based on a chiral 
Lagrangian that includes in addition to pions and nucleons, the A(1232) and the p and a 
mesons. The potential obtained from this Lagrangian consists of s- and -u-channel A and 
A pole diagrams as well as t-channel p and a exchanges. Convergence of all integrals is 
guaranteed by the use of cutoff fuctions associated with each vertex. Two different parame- 
terizations of the cutoff functions were considered: in model I the cutoff function was taken 
to be a product of form factors depending on the four-momentum squared of each particle 
present at the vertex. In model II the cutoff function was taken to depend only on the pion 
four-momentum squared. The parameters of the potential were adjusted to fit the empirical 
s- and p-wave phase shifts up to a pion laboratory kinetic energy of 360 MeV. Both models 
give good fits to the nN scattering data, and the resulting coupling constants are consistent 
with the commonly accepted values extracted from other observables. 
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While most of our results were for the Rarita-Schwinger A propagator, we compared the 
results of fits performed using the Rarita-Schwinger and Pascalutsa propagators for the case 
of the model I form factors. The differences in the coupling constants obtained suggests that 
a complete understanding of the baryon resonances with higher-spin is not possible without 
having unique higher-spin propagators. 

The good fits to the on-shell nN data for energies below 360 MeV suggests that a model 
of the ttN interaction based on the Bethe-Salpeter equation could form the basis for the 
analysis of pion photoproduction by the proper U(l) gauging of the Lagrangian |3~Ifl , and 
the analysis of the baryon resonances near and above the threshold for pion production. 
In fact, by extending the calculations to pion energies up to 600 MeV, we observe that 
in partial waves that don't have large inelasticities and don't exhibit evidence of baryon 
states, both our models give good representations of the data. For partial waves with large 
inelasticity there is also evidence for baryon resonances, and as a result we need to include 
first the coupling to inelastic channels, and then include bare baryon states into the potential 
if necessary. 

Although we have not included full three-body unitarity, which would require the dressing 
of the nucleon propagators and tcNN vertices, we have included those contributions to three- 
body unitarity resulting from the fact that we have not carried out any three-dimensional 
reduction, and as such our potentials depend on the relative energy. The inclusion of the 
dressed nucleon in the ttN intermediate states will increase the number of coupled channels, 
and is under investigation. 

By calculating off-mass- shell ttN amplitudes using the Bethe-Salpeter equation we can 
examine the low energy theorems, and study the questions associated with the analytic 



continuation of the physical tiN data to the Cheng-Dashen point [B4|, and the changes in 



the ttN sigma term as one goes from the Weinberg point, where the four-momenta of the 
pions is zero, to the Cheng-Dashen point, where the pions are on-mass-shell. This question 
is presently being examined and could shed some light on the inconsistencies between the 
"observed" £jv and cr/v as extracted from QCD models. 

It is clear from the work presented here that constraints on the coupling constants and 
form factors need to be improved before it can be established whether or not the present 
potential includes all the physics of ttN scattering at low energies. We are of the opinion 
that such constraints, particularly on the coupling constants, should come from QCD or 
QCD-based models. 
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cutoff mass minimum value 

A a rriN + 3m,r 

A w (mAr + 3m 7r )/2 

K p (m N + m w )/2 

K (m N + m 7r )/2 



TABLE I. The minimum allowed values of the cutoff masses in terms of the nucleon and pion 
masses. 



coupling constants 


RS 


Pas 




13.5 


13.5 




1.80 


12.1 


0.365 


0.741 




0.37 


0.193 


XA 


-0.11 






2.88 


2.73 




2.66 


4.11 


ScTTirdcNN /47T 


-0.41 


-3.80 


masses 


m N 


0.939 


0.939 


(0) 


1.34 


1.72 


mA 


1.232 


1.232 


(0) 


2.305 


2.60 




0.138 


0.138 


m p 


0.769 


0.769 


m a 


0.65 


0.69 


A N 


3.17 


4.90 


A A 


4.56 


3.20 


A^ 


1.77 


1.76 


Ap 


3.67 


3.06 


K 


1.30 


4.26 



TABLE II. The coupling constants and particle masses obtained using the Rarita-Schwinger 
(RS) and Pascalutsa (Pas) A propagators. The quantities in boldface were varied in the fits. All 
masses are in GeV. 
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l 2I2j 

Sn 
S31 
Pu 

Pl3 
P31 
P33 



BS (RS) 
0.177 
-0.101 
-0.083 
-0.032 
-0.041 
0.178 



BS (Pas) 
0.172 
-0.105 
-0.058 
-0.031 
-0.041 
0.187 



SM95 
0.175 
-0.087 
-0.068 
-0.022 
-0.039 
0.209 



KH80 
0.173 
-0.101 
-0.081 
-0.030 
-0.045 
0.214 



TABLE III. Scattering lengths and volumes obtained from the BS equation in units of m n 
compared to results from the SM95 [33] and KH80 [58] irN partial wave analyses. 



(2£+l) 









K p 


equation 


ref. 


0.365 


-0.11 


2.88 


2.66 


BS (RS) 


this work 


0.741 




2.73 


4.26 


BS (Pas) 


this work 


0.35 


-0.3 


2.5 


3.7 


ET 





0.43 




2.85 


1.8 


ET 


@ 


0.36 


-0.12 


3.13 


2.25 


Sm 


Q 


0.36 


-0.41 


2.90 


1.44 


BbS 





0.40 


-0.21 


3.36 


6.6 


Ka 





0.36 


-0.31 


3.03 


3.16 


Tr 





0.36 




3.1 


2.7 


Tr 






TABLE IV. Comparison between coupling constants obtained from the Bethe-Salpeter (BS), 
equal-time (ET), smooth (Sm), Blankenbecler-Sugar (BbS), Kadyshevsky (Ka) equations, and 
tree-level (Tr) calculations. 
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coupling constants 


n n = 2 


n-n = 4 


= 10 


model I (RS) 


91nn/^ 


13.5 


13.5 


13.5 


13.5 




4.23 


4.68 


5.98 


1.80 


0.365 


0.365 


0.371 


0.365 


f(°)2 / 4vr 


0.17 


0.20 


0.196 


0.37 




-0.13 


-0.24 


-0.18 


-0.11 




2.67 


2.63 


2.80 


2.88 


ftp 


2.18 


2.03 


2.15 


2.66 




0.86 


0.39 


0.48 


-0.41 


masses 


0) 

my 


1.18 


1.14 


1.11 


1.34 


(0) 


1.495 


1.492 


1.435 


2.305 


m<j 


0.88 


0.62 


0.64 


0.65 




1.34 


1.85 


2.73 


1.77 



TABLE V. The coupling constants and particle masses resulting from fits to the ttN data using 
different values of n n . The parameters from the model I fit using the Rarita-Schwinger propagator 
are also shown for comparison. The quantities in boldface were varied in the fits. All masses are 
in GeV. Particle masses not given are the same as those given in Table 0. 





n n = 2 


n-n = 4 


n n = 10 


A? 


0.874 


0.868 


0.822 


A^ 


2.47% 


2.51% 


2.79% 



TABLE VI. The renormalized pion cutoff masses A^ (in GeV) and values of A^ (expressed as 
percentages) . 
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(b) 





(c) 





+ 




(d) 



FIG. 1. Graphical representation of: (a) the BS equation for the full T-matrix, (b) the potential 
and (c) the T-matrix in terms of non-pole and pole parts, and (d) the BS equation for the non-pole 
T-matrix. 






FIG. 2. The equation for the dressed ttNB vertex. 





FIG. 3. The baryon self-energy. 
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FIG. 4. The diagrams included in the potential of the BS equation. 
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FIG. 5. The phase shifts obtained from the BS equation shown versus the pion laboratory energy, 

using the Rarita-Schwinger A propagator ( ) and Pascalutsa A propagator ( ). Data points 

from the VPI SM95 partial wave analysis are also shown. 
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FIG. 6. Contributions to the irN phase shifts as each diagram is added to the potential, in the 
following order: w-channel N pole ( — ), s-channel N pole (--), ^-channel a exchange (•••), t-channel 
p exchange (- •), u-channel A pole ( — ), and s-channel A pole ( ). 
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FIG. 7. The phase shifts obtained from the BS equation are shown versus the pion laboratory 
energy up to 600 MeV, for models I ( ) and II ( — ). 
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